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Removal and mixing of the coronal gas from satellites in 
galaxy groups: cooling the intragoup gas 
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ABSTRACT 

The existence of an extended hot gaseous corona surrounding clusters, groups and 
massive galaxies is well established by observational evidence and predicted by current 
theories of galaxy formation. When a small galaxy collides with a larger one, their 
coronae are the first to interact, producing disturbances that remove gas from the 
smaller system and settle it into the corona of the larger one. For a Milky- Way- 
size galaxy merging into a low-mass group, ram pressure stripping and the Kelvin- 
Helmholtz instability are the most relevant of these disturbances. We argue that the 
turbulence generated by the latter mixes the material of both coronae in the wake of 
the orbiting satellite creating a "warm phase" mixture with a cooling time a factor 
of several shorter than that of the ambient intragroup gas. We reach this conclusion 
using analytic estimates, as well as adiabatic and dissipative high resolution numerical 
simulations of a spherical corona subject to the ablation process of a constant velocity 
wind with uniform density and temperature. Although this is a preliminary analysis, 
our results are promising and we speculate that the mixture could potentially trigger 
in situ star formation and/or be accreted into the central galaxy as a cold gas flow 
resulting in a new mode of star formation in galaxy groups and clusters. 

Key words: hydrodynamics turbulence - cooling flows -galaxies: evolution - groups 
- interaction - methods: analytical - numerical 



1 INTRODUCTION 

The classical theory of galaxy formation establishes that 
gas follows dark matter in its non-linear evolution as it 
forms virialized dark matter haloes. In this process, the 
gas is heated by shocks and adiabatic compression, acquir- 
ing the virial tempe rature of the halo it has fallen into 
I White fc Rees|[l97Sl '). Afterwards, the gas can lose energy 
effectively through radiative cooling if the cooling time is 
much smaller than the free-fall time (t coo i < iff), produc- 
ing a contraction of the gas towards the centre: a cooling 
flow that fuels star formation. The critical mass for effective 
cooling depends on the density, temperature and metallicity 
of the gas. Assuming constant gas density and temperature, 
one can show that the gas in haloes with primordial gas 
masses > 10 n M o cannot cool effectively (e.g. see Fig. 8.6 of 
iMo et alJlioich ; this is why galaxy groups and clusters are 
expected to have a corona of hot gas at the present time. 

X-ray observations of relaxed clusters and groups show 
indeed the presence of an extended gas corona within the 
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intra-cluster and intra-group medium (hereafter called ICM 
in both cases), which is approximately in hydrostatic equilib- 
rium within the gravitational potential of the cluster/group. 
For the most massive systems, the density profile of this dif- 
fuse gas approximately t races that of the pre dicted dark 
matter distribution (e.g. IVikhlinin" et alj |2006h Since the 
dominant cooling channel in these systems is bremsstrahlung 
radiation, then t coo i oc T 1,/2 p _1 (for primordial composi- 
tion) which implies that the cooling time scales are typically 
shorter in the central regions, and hence, in the absence of 
a balancing heat source, a cooling flow would develop in the 
centre. Close to the virial radius of the system however, the 
typical cooling times are larger than the Hubble time. 

A hot extended corona is also predicted to exist within 
the haloes of massive gal axies, including those at the scale 
of the Milky- Way (MW; I White fc Frenkl Il99lh . In the ab- 
sence of additional energy sources, and under idealized cir- 
cumstances, the majority of this gas would be expected to 
cool and collapse into the halo centre over a Hubble time; 
in contrast, roughly 20% of the expected baryons in these 
haloes are o bserved to be in a cold phase (stars and cold 
gas, see e.g. iFukugita fc Peebledliooi . l2006t and references 
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therein). This overcooling problem can be avoided by in- 
voking feedback processes that prevent the gas from cooling 
catastrophically. Recent cosmological simulations that in- 
clude radiative cooling, star formation and feedback predict 
that MW-size haloes should have a gas co rona surviving to- 
day in quasi-hydrostatic equi librium (e.g. ICrain et al.ll201ol ; 
Ivan de Voort fc Schave|[2012l '). It has also been argued that 
even without feedback, the combined effect of buoyancy and 
thermal conduction suppresses the thermal instability, and 
thus the formatio n of cold clouds, in stratified co ronae in 
MW-size galaxies (|Binnev et al.ll2009l ; lNipotill2010t ). 

The presence of extended coronae has been confirmed 
for elliptical galaxies through X -ray observations (see review 
by iMathews fc Brighentll l2003h but the evidence for MW- 
like galaxies is not conclusive. Recent observations have de- 
tected diffuse X-ray emission around disc-dominated galax- 
ies, but is approximately two orde rs of magnitude lower 
than the e xpectations ba sed on the IWhite fc Frenkl (| 199 if ) 
model (e.g. [Li et al1l2007r ); supernovae he ating is advocated 



as the likely cause of this discrepancy (e.g. Crain et al ] |20ld : 



lAnderson fc Bregmarj|201ll ; iDai et "alj|2012h . We note that 
although there is evidence for the presence of gaseous coro- 
nae around galaxies, their cosmological origin as a shock- 
heated remnant is not settled. There are other possible 
origins such as enriched gas produced by galactic stellar 
feedb ack (e.g. lCiotti fc Ostrikerlll997l ; lMathews fc Brighentil 
Il998h . 

It is thus relevant to consider in detail the impact of the 
gaseous coronae on the process of galaxy formation, partic- 
ularly during galaxy mergers. In these events, the gas in 
the merging galaxies is subject to different processes that 
remove gas from the satellite: i) ram p ressure stripping re - 
moves the most loosely bound material l|Gunn fc Gottlll972[ ) 
creating a wake that trails behind the satellite (recent obser- 
vations sho w evidence of suc h wakes in galaxies falling into 
clusters, e.g. lSun et al.|[201ol . but the stripped gas in the ob- 
served cases is the cold ISM of the galaxy, not its corona); ii) 
thermal cond uction causes evaporat ion of the colder gas in 
the satellite (|Cowie fc McKedll977l ); hi) if the flow around 
the boundary interface is la minar then v iscosity will cause a 
drag that will remove gas (|Nulsenlll982l ); iv) fluid instabil- 
ities such as the Kelvin-Helmholtz (KH) instability caused 
by the velocity shear between the layer interface also lead 
to gas removal. 

All these processes have been studied in the past and for 
the case of galaxies merging into clusters, particular atten- 
tion has been given to ra m pressure stripping using analytic 
and numerical methods IIMori fc Burke rt 2000; Heinz et al.l 
120031 ; iMcCarthv et al.l 120081 '). Most of these studies focus 
on the impact of the removal of gas in the star formation 
efficiency of the satellite: by completely stripping the gas 
corona, a fuel reservoir is removed, resulting in a quenching 
of star formation relative to a case where the galaxy evolves 
in isolation (this process is usually called "strangulation"). 
Recent analyses show that a substantial fraction of the 
corona can actually survive if the satellite is massive enough, 
and through su bsequent radiative coo l ing, can contribute to 
star formation (|McCarthv et al.ll2008l : iMoster et al.ll201lT ). 

The fate of the stripped gas has not been analysed in 
detail however. As mentioned before, the gas from the satel- 
lite removed by ram pressure will form an elongated wake 
trailing the satellite. Due to the gravity of the satellite, a 



similar wake will be formed with gas from the ICM that 
gets gravitationally focused l|Stevens et al.l 1999: S akellioul 
2000). Thus, in general we would expect a wake composed 
of material from both the ICM and the corona of the satel- 
lite, implying that both gaseous components mix during the 
removal processes. This mixing is particularly efficient dur- 
ing the turbulence caused by KH instabilities. If the gas 
removal rate is large and the mixing process is efficient then 
the mixture will have a temperature that might be signifi- 
cantly lower than the temperature of the ICM. This would 
imply that the cooling time in the wake is lower than that 
in the ambient medium, making the mixture a potential fuel 
source for star formation. This i dea has been proposed re- 
cently by iMarinacci et alj (|2010t ) in the context of a mode 
of gas accretion triggered by the interaction of Hi clouds 
expelled from galactic discs in a galactic fountain, and the 
corona of the galaxy (see also lHeitsch fc Putmanl [20091 ). In 
this paper we argue that the same mechanism is at play 
when the coronae of the merging galaxies interact. 

Most of the hydrodynamical simulations of galaxy for- 
mation carried out to date have used the smoothe d particle 
hydrod ynamics (SPH) technique (for a review see ISpringell 
l2010bl ). This is a formulation that has excellent conservation 
properties in terms of mass, energy, momentum and entropy 
and it is fully Galilean invariant. However, SPH consider- 
ably suppresses the formation of turbulence triggered by in- 
stabil ities and thus it i s not able to account fo r gas mix- 
ing l|Agertz et al.ll2007l ; [Bauer fc Springeill2012r ). Adaptive 
Mesh Refinement (AMR) codes based on an Eulerian formu- 
lation are the other method currently used for simulations of 
galaxy formation. Although they suffer their own shortcom- 
ings (like the lack of Galilean invariance), perhaps the most 
important difference with SPH methods is precisely that 
of th e treatment of turbulence and mixing (|Mitchell et al.l 
[2009j), which are not suppressed in AMR. 

In this paper we take a first step in investigating the 
effect of mixing in the interaction of coronae in galaxy merg- 
ers, focusing on the formation of the trailing wake. For this 
purpose we start by revisiting analytically the most impor- 
tant gas removal processes and estimate the properties of the 
resultant gas mixture (Section [3J. We then verify these esti- 
mates using 3D adiabatic (non-radiative) , as well as dissipa- 
tive simulations of the motion of a spherical corona through 
a constant density ICM using the AMR code FLASH (Sec- 
tion Our results are presented in Section [4] and finally a 
discussion and our conclusions are given in Sections [5] and 
[S] Although our simulations lack several important physical 
processes that could impact the results, our objective here is 
to make an initial assessment of whether or not the removed 
coronal gas could continue to be a source of star formation. 



2 THEORETICAL EXPECTATIONS 

In this section, we present analytic estimates of what might 
happen to the gas in the corona of a satellite as it merges 
with that of a group/cluster. The corona is embedded in a 
spherical dark matter halo which is assumed to have a ra- 
dial density distribution given by the Navarro-Frenk- White 
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Table 1. Typical values for the global properties of a MW-size galaxy, a 
low-mass group and a cluster. 



System Af v (M Q ) r v (kpc) V v (km/s) T V (K) c s (km/s ^ 



MW-size galaxy 


10 u 


200 


146 


7.7 x 10 5 


135 


low-mass group 


10 13 


430 


315 


3.6 x 10 6 


291 


cluster 


10 14 


931 


680 


1.7 x 10 7 


627 



a The local sound speed is given by c 2 = yP/ p with 7 = 5/3. 



(NFW) profile |Navarro et al .11 19961 . [19971 1*1 which in terms 
of the dimensionless variable s — r/r v , where r v is the virial 
radius, can be written as: 



p(s) = 



Pcrit 5c 



(1) 



cs(l + cs) 2 ' 

where c is the concentration parameter, p cr it is the critical 
density and S c is a function of the concentration: 

_A 

c ~ 3"log(l + c)-c/(l + c)' () 

The virial radius r v is defined as the radius where the mean 
density of the halo is equal to A times the critical density: 
M v /(47rrv/3) = Ap cr i t , where M v is the total mass within 
r v . The value of the density contrast A is conventional and 
varies in the literature, but we use A = 200 for this work. 
The concentr ation correlates strongly with virial mass and 
redshift (e.g. iGao et aill2008l ). When a numerical value is 
needed, we adopt the mass-concentration relation found in 
iNeto et al.l (|2007t ) b ased on the Millenn ium simulation of 
structure formation (|Springel et al.ll2005h . 

We further assume that the dark matter haloes con- 
tain a distribution of gas that also follows the NFW profile 
but normalised by the universal fraction of baryons to dark 
matter /b =■ Qb/{^m — Qb). Whenever it is needed, we take 
fi m = 0.25 and fib = 0.045 as the values for the contribu- 
tions from matter and baryons to the total energy density 
of the Universe, respectively. 

The pressure profile is obtained by assuming that the 
coronal gas is in hydrostatic equilibrium: 

dP d<t> 

where </>(s) is the total gravitational potential (dark mat- 
ter + gas). To obtain a boundary condition for Eq. @, 
we note that the virial velocity of the dark matter halo 
Vv = GM v /r v determines a characteristic temperature for 
the corona T v = /imirK, 2 / (2&b), where fi is the mean molec- 
ular weight of the gas (/i _1 = 1.71 for gas with primordial 
composition), run is the mass of a hydrogen atom and ks is 
the Boltzmann constant. This characteristic value, and also 
the assumption of hydrostatic equilibrium, are motivated by 
the formation of the corona when the gas fell into the halo 
and got shock- heated to temperatures of the order of T v . The 
boundary condition for Eq. ([3]) is then set by P v given by the 



1 The most recent N-body simulations of structure formation 
suggest that the density profile is actually shallower than NFW 
towards the centre with an Einasto profile giving a better fit 
(Springcl et al. 2008). However, since we are not interested in the 
central regions, the NFW profile is an appropriate choice. 



value of the virial temperature of the halo and the equation 
of state (EoS), assumed to be that of an ideal gas. Since T v 
is actually close to the peak of the temperature profile, we 
apply a small correction (a factor of ~ 3/4) to the value of 
P v to obtain a temperature profile that decays as a power 
law at large radii, which is a behaviour seen in observations 
of ne arby galaxy clusters and groups (e.g. IVikhlinin et al.l 
2005), as wel l as in full hydrodyn amical simulations of such 
systems (e.g. iBorgani et ahll2004l '). 

Unless otherwise stated, we assume that the satellite 
halo enters the ICM with a relative velocity u sa t which is of 
the order of the virial velocity of the host V^host- Through- 
out this section, we compare the relevance of different pro- 
cesses during the interaction of the hot coronae in two cases: 
a MW-size galaxy falling into a low-mass group (galaxy- 
group collision) and the same galaxy falling into a cluster 
(galaxy-cluster collision) . The merger mass ratio in the cases 
of interest is ^ 10, so the gravitational potential of the most 
massive system is not significantly affected. The character- 
istic parameters of these three systems appear in Table [1] 



2.1 Removal of gas from the satellite 

2.1.1 Ram Pressure Stripping (RPS) 

As the satellite moves through the ICM, the ram pressure 
(fram = Phostfsat) imparted to its corona will be sufficient 
to remove gas from the satellite if it can overcome its gravi- 
tational restoring force per un it area l|McCarthv et al ] |2008l . 
based on lGunn fc Gott|]l972l ): 



> a. 



Gm to t(< r)p sat (r) 



(4) 



where p sa t(r) is the gas density profile of the corona of the 
satellite, m to t(< r) its total mass profile (dark matter + 
gas), and a g is a geometric factor of 0(1) depending on 
the total mass and gas density profiles of the satellite; for 
instance, a g = 7r/2 for a singular isothermal sphere. The 
radius where the lhs and rhs of Eq. (Q are equal is the pre- 
dicted radius r rps that the gas sphere will have after strip- 
ping; m sat (< r rps ) is then the surviving mass. 

As long as the relative velocity between both sys- 
tems is supersonic, the collision drives a shock that prop- 
agates through the satellite; otherwise a compression wave 
is formed that propagates with the local sound speed. It is 
possible to estimate the speed of this shock and its subse- 
quent propagation as a function of the radius of the satellite 
v s h(r) by solving a one-dimensional shock tube problem, ne- 
glecting the effects of gravity (see Appendix [XJ. The char- 
acteristic RPS time scale is then roughly given by the time 
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it takes for this shock to cross the diameter of the satellite 

dr 



trps ~ 2 . 

Jo v M r ) 



(5) 



The average mass removal rate from RPS is then given by 

frps) ^satf*^ ?"v,sat) 



(6) 



2.1.2 Tidal Stripping (TS) 



Tides remove gas and dark matter from the satellite as it 
falls into the host. The tidal radius r t , defined as the ra- 
dius where the external differential tidal force from the host 
exceeds the b inding force of the sa tellite, is approximately 
given by (e.g. I Taylor fc Babulllaxnl ): 

3 _ Gm to t(< r t ) 



11 w^ + G[2M{<r)/r^-4-K Phoat (r)]' (7) 

where w is the angular speed of the satellite located at a 
distance r from the host, which has a mass and density pro- 
file M(< r) and phost(r) respectively. The material outside 
r\ is tidally removed from the system within a time scale 
its = torb/A where t or b = 2n/w is the instantaneous orbital 
time and A is an efficiency parameter. It is important to 
note that Eq. (0 is only a rough approximation^. The un- 
certainties associated with this equation are parameterised 
by the factor A that has values between 1 and 6 according 
to numerical si mulations of dark matter subhaloes orbit- 
ing host haloes ilTavlor fc Babulll200ll ; IZentner et ai1l2005l ; 
iDiemand et al. 2007 ). We ta ke a fiducial value of A = 3.5 
as found bv lGan et ail l|2010h , and define the mass removal 
rate due to TS as: 



rn tB 



m sat (< r t ) — m sat (< r v , sa t) 
tZ ' 



(8) 



Besides RPS and TS, the corona of the satellite can be 
subjected to other mass loss m echanisms wh ich we mention 
subsequently (closely following lNulsen l ll982l ) . 



2.1.3 Laminar viscous stripping 

If viscosity is relevant in the ICM, then there will be a drag 
in the corona of the satellite resulting in a gas mass loss 
which is roughly given by (ignoring gravity): 



12 \ 2 

) 7IT 8a t/9host«sat, 



(9) 



where r sat is the instantaneous maximum radial extent of 
the corona of the satellite and Re is the Reynolds number: 



Re = 2.8 



(10) 



v Ahost J VC^host^ 

with Ahost being the effective mean free path of ions in the 
ICM. In the absence of magnetic field s and for Coulomb 
collisions this is given by (|Spitzerl ll956): 

2 / ; \ -1 



Ahoat W 11 kpC 



Phost/rriH 



(11) 



10 8 K/ V10" 3 cm- 3 
For Eq. (O to be applicable, the gas flow needs to be 

2 For instance, the tidal radius is only well-defined for circular 
orbits, whereas realistic orbits are in general eccentric. 



laminar around the boundary of the satellite, i.e., Re ^ 30, 
and also A host < r sat . Since T^ ost /p host oc r 1-9 for the NFW 
profile (and for the temperature given by the EoS assum- 
ing hydrostatic equilibrium), Ahost (r) < Ahost (rv.sat)- At the 
virial radius of the host, the mean free path is of 0(1 kpc) 
and since « sa t/ c s,host > 1 (the initial shock is transonic), 
the Reynolds number is then very large and viscosity is ir- 
relevant (see Table [2j). Due to RPS and TS, r sat will be 
reduced as the satellite spirals inwards to the centre of the 
host but the ratio r sa t /Ahost will certainly remain very large. 
The changes in the ratio Vs a t/c s ,host will depend strongly on 
the orbit but will likely remain larger than 1, increasing to- 
wards pericentre, and dropping below 1 near apocentre. If 
apocentre is near the virial radius of the group (cluster) and 
the radius of the satellite has been reduced by a significant 
factor, then Re could drop below 30 in the cluster environ- 
ment. This is unlikely to happen in the group environment 
because the mean free path of the ions in the ICM is an 
order of magnitude smaller than in the cluster environment. 
We neglect viscous stripping in the remainder of this work. 



2.1.4 Evaporation due to thermal conduction 

The rate of evaporated mass due to thermal conduction from 
a spherical cloud embedded in a hot tenuous gas, neglecting 
gravity, radiation, ionization and magnet ic fields, is approx- 
imately given by l|Cowie fc McKedl 19771) : 



47T7"g at phost Cs.host <l>sF(a ), 



(12) 



where <j) s is a parameter of 0(1) and 00 = 1. 84 Ahost /fsat^s 
determines when the regime of classical conductivity still 
applies. This happens when the mean free path is short com- 
pared to the temperature scale height Lt = T/\ VT| between 
the boundary of the cloud and the hot gas surrounding it 
(Ahost <C Lt); in this regime, F(ffo) = 2ao- For the set of in- 
teracting systems we are considering, r sa t /Ahost > 0(100), 
thus, 00 < 0(0.01), which makes the mass loss rate due 
to thermal conduction much lower than that of the other 
processes for the group environment (see Table [2} . For the 
cluster environment, evaporation is still subdominant but it 
cannot be entirely neglected. 

2.1.5 Turbulent Stripping (Kelvin-Helmholtz Instability) 

If the inertial forces in the surface layers between the two 
gases are larger than those related to viscosity (i.e. if Re > 
30), then the velocity shear between the layers will drive 
KH instabilities. The collision between the coronae of the 
satellite and the ICM also drives a reverse shock that prop- 
agates into the ICM (contrary to the forward shock prop- 
agating into the satellite and related to RPS). Since, typ- 
ically, Wsat ~ Cs,host at infall, this shock is generally weak. 
The post-shock flow in the vicinity of the boundary of the 
satellite, where the instabilities develop, will then be sub- 
sonic (propagating at a speed v ps , see Eq. A2). Under this 
assumption, and taking into account the compressibility of 
the gas flow, the mass loss rate for complete stripping due 
to KH instabilities is: 



771KH ~ Kr sat p host V ps . 



(13) 



In this equation, the dominant wavelength of the pertur- 
bations is of the order of the radius of the corona of the 
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Table 2. Typical values for the parameters defining different processes of gas removal for a MW-size galaxy falling 
into a group/cluster (see Table[TJ. Times are given in Gyr and mass loss rates in Mo/yr. 



Environment 


^rps/^v 


^rps 




n/r v 


Us 


mta 


^host/ r v 


Re 


*kh 


rn K u 


T^cvap 


galaxy-group 


0.78 


2.1 


15 


0.65 


2.4 


22 


7 x 10~ 4 


4800 


1.2 


39 


0.5 


galaxy-cluster 


0.45 


1.9 


50 


0.65 


2.4 


22 


1.4 x lO" 2 


250 


0.56 


53 


23 



satellite. This is because the stripping process begins with 
perturbations that have shorter wavelengths; as they de- 
velop, they smooth the velocity gradient between the layers 
making the interface thicker. This damps modes with wave- 
lengths smaller than the thickness of the interface, while 
larger wavelength modes continue growing, stripping more 
material and making the interface thicker. This process con- 
tinues until the layer has a thickness of order r sat . 

As we mentioned before, the velocity flow near the layer 
is expected to be subsonic. We estimate the values of v ps by 
solving a one dimensional shock tube problem (see Appendix 
|A"{ using the value of the hydrodynamical variables at the be- 
ginning of the interaction (assuming initial pressure equilib- 
rium between both systems). We find v ps ~ 208 (260)km/s 
for the group (cluster) environments resulting in mass loss 
rates which are of the order of those by RPS and TS (see Ta- 
ble [2j . The time scale for the formation of the largest wave- 
length instabilities can be roughly estimated for an ideal case 
where the fluids are incompressible, the surfaces are fiat, 
gravity is neglected, and the overdensity near the boundary 
is in t he linear regime: S = p sa t/phost — 1 <C 1. In this case 
(e.g., iDrazin fc Reidlll98ll ): 

*KH~— , (14) 

Vsat 

which is of the order of the sound crossing time. 

Gravity and density stratification of the gas are stabi- 
lizing effects acting against the development of KH insta- 
bilities. By assuming that the interface has a sharp density 
discontinuit y between the gas in the satellite pi and the 
ICM gas pi, iMurrav et all (jl993h found that the dominant 
wavelength for gas ablation, k — 27r/r sat , is stable if: 



9 > 



27vp 1 p 2 v S: 



r sa t (pi 



(15) 



where g = Gm to t(< r sa t) /r^tA i s t ne gravitational acceler- 
ation at the interface. This condition can be rewritten in 
terms of the gravitational restoring pressure (Eq.[3]): 



restore \< s a t 



(r B . 



> 



27ra g/ 0iu 5at 
1 - {pilpif 



(16) 



The rhs of Eq. I|16l) is expected to be much larger than the 
restoring force in the outskirts of the satellite, and thus, 
the corona will not be stabilized by its own self-gravity but 
the mass loss rate in Eq. (|13[) is probably going to be re- 
duced. Other physic al processes, such as magnet ic fields (e.g. 



Malag oli et al.lll996t ) and radiative cooling (e.g. IVietri et al.l 
1997h might stabilize the corona of the satellite. 

Given the parameters of the interacting galaxies, it 
seems likely that RPS, TS and KH instabilities will act on 
similar time scales, competing to remove the gas from the 
corona of the satellite. Besides gas removal, all these pro- 
cesses will also mix the stripped gas with the gas in the 
ICM with important consequences as we discuss below. 



2.2 Cooling due to mixing 

As the satellite orbits the host, it los es energy and angular 
momentum due to dynamical friction l|Chandrasekharlll943r ) 
eventually merging completely with the host. The typical 
time scale for this to happen is ~ 4(23) Gyr for the satel- 
lite merging into the group (cluster) (following the study of 
iBovlan-Kolchin et al.l 120081 ). This implies that the satellite 
will remain in orbit for a reasonable time before getting too 
strongly disrupted. Thus, the gas removed from the satel- 
lite will have enough time to mix with the ICM, changing 
its properties. In the following, we consider this gas mix ing 
process, inspired by the study of iMarinacci et alJ ( |2010h . 

The coronal gas is being removed at a rate given by the 
different stripping processes. If we take the mass loss rate 
by RPS as a reference, then m Icm = a IC mm lps , with a rcm > 
0(1) representing the mass loss for other processes. The 
removed gas mixes with the ICM (possibly quite efficiently 
due to the turbulent flow generated by KH instabilities) and 
falls behind the satellite forming a wake. In the time it takes 
the satellite to move a distance s, the gas mixture occupies 
a volume V w = A v s, where A w is the cross section of the 
wake. If we neglect the compressibility of the gas, then the 
mass from the corona and the ICM contained within V w 
is: M w = Mhost, w + Msat.w = Avsphost + m rcnl s/v sa t. If the 
gases in the corona and the ICM have the same composition, 
then the temperature of the mixture is just: 

Afhost.wThost + Af S at,wT S at _ 



-^host.w + -Msat,w 
j4wphost.Thc.st ~\~ ?i^rem?s a t / Vsa.t 
^4-wphost ~~T" Tflrem / V B at 



(17) 



Note that thermal conduction is necessary to reach this tem- 
perature, but this can happen efficiently on small scales due 
to turbulent mixing. Depending on the relative difference 
between the gas from the satellite being deposited in the 
wake fhrem/fsat to that of the ICM j4 w phost, the temper- 
ature of the mixture can decrease significantly relative to 
Tkost, down to a minimum of T sat . 

The characteristic scale for the cross section of the wake 
is given by a combination of the gravity of the satellite, 
the inertia of the stripped gas and turbulent effects. For 
simplicity, we parametrize it as A w — fi-Kr'i, where r B is 
the scale radius of the halo of the satellite: r s = r VjS at/c 
(with c being the concentration of the halo), and ft is a free 
parameter which in the absence of turbulence is mainly set 
by the velocity of the satellite^. We note that in reality, the 



3 This cross section seems to be a good guess according to the 
results we found for a low resolution simulation analogous to the 
non-radiative simulation described in Section[3]for a subsonic col- 
lision (u S at = 150km/s). We found that the cross-sectional radius 
of the wake separating the gas in the ICM that flows without 
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Figure 1. Cooling time as a function of radius for a group/cluster 
with an NFW density profile under hydrostatic equilibrium. The 
thick solid lines are for an isolated system; the thin solid lines 
show what would happen if gas in the ICM at different radii 
mixed with coronal gas removed from an infalling satellite (T sat = 
7.7 X 10 5 K) moving with a constant speed equal to the virial 
velocity of the group/cluster (the cooling time refers in this case 
to the gas in the mixture). The gas is removed at a constant 
rate m rcm = a rcul rh TpB and mixes with the ICM in a region with 
a cross sectional area A w = /3irr B . The thin solid lines are for 
c*rcm = /3 = 1. The dashed (dotted) lines are for c* rcm = 1/2(2) 
and the long-dashed (dotted-dashed) lines are for /3 = 10(1/10). 



gas mixes across layers of different sizes in the turbulent flow 
creating a broad distribution of j4 w values. 

The time scale for gas to lose its thermal energy due to 
radiative losses is given by: 



icool 



3/2nk B T 
n e rm Accci (T, Z) ' 



(18) 



where n = phost/(M m H), n a and nrr are the number densities 
for electrons (n c = 0.52n) and hydrogen, respectively; and 
Acooi (T, Z) is the cooling function that depends on temper- 
ature and metallici ty. We obtain the cooling function from 
the tables given bv lGnat fc Sternberg! (|2007l ). 

Fig. |l} shows the cooling time for the gas mixture as 
a function of the distance to the centre of the group (clus- 
ter) assuming that the ICM in the host mixes with gas from 
the satellite with a single temperature (T S at = Tv^at). We 
also assume that the density of the mixture is the same 
as the density of the hoslQ, and that the mass loss rate 
and velocity of the satellite remain constant with radii. We 
have used a metallicity of \ogio(Z / Zq) = —0.5, which is 
representative of the values found in groups and clusters 



mixing with the one that mixes is ~ r B , i.e., = 1. For higher 
velocities, /3 is likely to increase asymptotically to (r v /r s ) 2 . 
4 Note that it is perhaps more reasonable to assume that the 
mixture will have the same pressure as the ambient ICM; in this 
case, its density will be higher and thus its cooling time will be 
lower than the one we have estimated. 



l|Rasmussen fc Ponmanl [2009). The thick solid lines show 
the reference cases where there is no mixing, in blue and 
red for the group and cluster, respectively. The thin blue 
and red lines show the case when mixing is effective, with 
Qrcm = /3 — 1 for the group and cluster, respectively. Only 
for the group, we show lines where the ratio rh rcm /i7 sat is 
decreased (increased) by a factor of 2, dashed (dotted) blue, 
and for changes in the parameter fJ increasing (decreasing) 
it by a factor of 10, long-dashed (dotted-dashed) blue. The 
former of these explores variations in the efficiency of gas 
removal, and the latter on the effective size of the region 
where material mixes. 

The regions where the different assumptions made to 
arrive at Fig. (flj are likely to be true, are those close to 
the outskirts of the host. We see that for r/r v> h os t ~ 0.6 the 
cooling time of the mixed material is of 0(1 Gyr) for the 
group environment, which is an order of magnitude lower 
than the cooling time of the ambient ICM. For the cluster 
environment, the reduction in the cooling time is also quite 
substantial but certainly above a few Gyrs. 

If a volume element of gas can cool down quickly enough 
by radiating its thermal energy, then, in the absence of any 
source of heating, it will lose all its pressure support and 
collapse in a free-fall time: 



iff = 



3tt 



32Gpiocal 



(19) 



The condition f coo i = tft separates the gas clouds that can 
cool down effectively from those that cannot. The local free- 
fall time for gas at r ~ 0.6r- Vj host in the group(cluster) envi- 
ronment is ~ 4Gyr. 

Therefore, from the theoretical expectations described 
in this section we conclude that: the coronal gas removed 
from the infalling satellite will mix with the ICM of the host, 
creating a mixture with a cooling time potentially lower than 
the local free-fall time. The net effect is an increase in the 
cooling efficiency of the gas mixture trailing behind the merg- 
ing satellite, relative to the ambient ICM. Of course, the gas 
mixture will have a cooling time which is larger than that 
of the initial corona. 

Note that we have chosen a mass for the MW-size 
galaxy which is at the low end of the range of Milky- Way 
halo masses allowed by current e stimates (see e.g. section 
5.1 of iBovlan-Kolchin et ai]|2012T ). which results in a virial 
temperature for the cor ona that is relatively low c ompared 
to other estimates (e.g. iFukugita fc Pe ebles 2006) estimate 
T v = 1.4 x 10 6 K since they assume M v = 2.6 x 1O 12 M ). 
A hotter corona would result in a reduction of the cool- 
ing efficiency of the mixed ICM gas. We have also chosen 
a low-mass group (10 13 M©), corresponding to a galaxy- 
group merger mass ratio of 1/10. Mergers with lower mass 
ratios ar e typically more comm on in a cosmological setting 
(see e.g. iFakhouri et alJ feoiO). and as implied by Fig. |T} 
for the case of the galaxy-cluster collision, a merger with 
a group of higher mass (virial temperature) would increase 
the temperature of the gas mixture. Nevertheless, the 1/10 
galaxy-group merger we are taking as a reference should be 
common enough to be of interest. We leave the analysis in a 
full cosmological setting for a future work and concentrate 
in the rest of this paper on the formation and properties of 
the wakes produced by the galaxy-group interaction. 

In the previous analysis, we neglected the compression 
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of the gas during mixing, which could in principle reduce 
the cooling time and the free fall time even further. How- 
ever, since an adiabatic compression would also increase the 
temperature of the gas (heating the mixture instead of cool- 
ing it), it might be that the two effects cancel each other. 
We answer this question in the following section. 

It is in principle possible that a fraction of the gas trail- 
ing behind the satellite is gravitationally re-accreted. How- 
ever, we find that for the parameters in Table[T] the presence 
of the satellite in the environment of the group(cluster) re- 
sults in a Bondi radius < 53(ll)kpc, which is smaller than 
the virial radius of the satellite. Since RPS happens fairly 
rapidly, the effective radius of the corona is quickly reduced 
while the mass of the dark matter satellite remains fairly the 
same at first since TS is expected to act more slowly (see 
Table[2J. Nevertheless, the stripping of the gas is not enough 
to reduce r sat to the Bondi radius. For this to happen in the 
group (cluster) medium, the stripping efficiency would need 
to be much higher, followed by a period of relatively low 
orbital velocities (low compared to the virial velocity of the 
host). The satellite could then potentially accrete material 
from the ICM. Although this situation is not common for 
the systems we are studying, in the case of the galaxy-group 
collision, a "fine-tuned" orbit could potentially have a Bondi 
radius smaller than the extent of the corona of the satellite 
near apocentre. 



3 SIMULATIONS 

We use the Eulerian code FLASH jFrvxell et aUl200Ch . that 
uses an AMR technique to solve the hydrodynamical equa- 
tions for the fluid. With this technique, FLASH is able to 
increase the resolution of the grid solver in the regions of 
most interest while keeping a coarser grid elsewhere. FLASH 
also includes the possibility of adding particles to the simu- 
lation in two varieties: active and passive. The former have 
mass and participate in the dynamics of the system, and we 
use them to simulate the dark matter halo. Passive particles 
have no mass and do not contribute to the dynamics, they 
are Lagrangian tracers of the flow and are useful to study 
the degree of mixing. The code uses a particle-mesh method 
to map the quantities from the mesh into the particles and 
vice versa. The Eu ler equations are solved usin g a piecewise- 
parabolic method (|Colella k, Woodwardll984h and the grav- 
itational pote ntial is solved using a parallel multigrid solver 
|Rickerll200sh . 



3.1 Initial conditions and setup 

The satellite consists of a spherical distribution of collision- 
less dark matter and gas with an ideal equation of state. 
Both of these fluids have a density profile given by Eq. (0. 
Dark matter is simulated with a set of N particles with 
their positions and velocities given by an algor i thm b ased 
on the method presented in Kaz antzidis et all (|2004l ). To 
assign the position of each particle, we randomly draw a 
value of the enclosed mass M(< r) in the interval [0, M v ] 
and then invert the function M(< r) to find the distance of 
that particle to the centre of the halo. The velocities of the 
dark matter particles are assigned from their actual energy 



distribution function give n by the Eddington formula (e.g. 
iBinnev fc Tremaindfl987l ): 

i r r d 2 p i ( do\ l 

where ^ — —(f) is the relative gravitational potential and e 
is the relative energy e = ^ — v 2 /2. We construct a table 
with values of e and /(e) that allows us to linearly interpo- 
late for any value of the relative energy. The magnitude of 
the velocity of each particle is then drawn at random from 
the distribution functi on 2 \^(r) — e] 1 ^ 2 /(e ) using a rejec- 
tion sampling method (|von Neumann! H 95 1). The direction 
of the position and velocity vectors is chosen at random from 
the unit sphere. This method ensures that the dark matter 
halo is a system in equilibrium (except of course from im- 
balances caused by discreteness effects). 

Gas is simulated by a hierarchy of cells using the AMR 
technique. The initial pressure in the cells is given by the 
condition of hydrostatic equilibrium as described at the be- 
ginning of Section[2] The velocities of the particles are set to 
zero. Finally the gas is set to have a mean molecular weight 
of fi' 1 = 1.71. 

In this way we generate an initial equilibrium configu- 
ration for a MW-size galaxy having the parameters given in 
Table [T] We have tested that this system remains in equilib- 
rium for at least ~ 5Gyr when evolved in isolation. Nonethe- 
less, two clear features appear in the density profile during 
this time: a small core (~ 5kpc) that develops in the centre 
due to the limited resolution of the grid, and a small fraction 
of the dark matter particles (~ 2%) expand to larger radii 
due to the tapering of the density profile at the virial radius. 

Finally, to simulate the ICM, the system is surrounded 
by gas with constant density (lCP 28 g cm~ 3 ) and temper- 
ature (3.6 x 10 6 K) moving relative to the satellite with a 
speed i>sat = K,, group = 315 km/s. The simulation box has a 
length of 1.9 Mpc with periodic boundary condition^]. The 
resolution of the simulation uses 32 cells for the coarsest level 
with 6 levels of refinement to give a maximum spatial reso- 
lution of ~ 1 kpc. We use 50 3 active particles to represent 
the dark matter halo of the MW-size galaxy, and 100 3 (50 3 ) 
tracer particles to sample the initial gas density distribution 
of the ICM (satellite). 

We consider two simulations, one without radiative 
losses and one where the gas is allowed to cool radiatively. 
This allow us to understand the impact of cooling in the for- 
mation of the turbulent wake: efficient cooling in the corona 
of the satellite prior to the removal of the gas would increase 
its density and reduce its sound speed, and hence suppress 
the KH instabilities, greatly reducing the amount of stripped 
gas. Central feedback is expected to reduce the cooling rate 
in the satellite, preventing a catastrophic collapse that would 
otherwise occur both in nature and in the simulations (this 
collapse also leads to prohibitive simulation time steps). To 
mimic the effects of feedback we multiply the cooling rate 
by a function of the distance to the center of the satellite: 
1—1/(1 + (r/r c ) 3 ) where r c — 0.15 r v . In addition, we in- 
troduce a temperature floor of 10 5 K, below which the gas 

5 Periodic boundary conditions are necessary to use the parallel 
multigrid solver for the Poisson equation in the version of FLASH 
we used. 



© 0000 RAS, MNRAS 000, 000-000 



8 Jesus Zavala, Michael L. Balogh, Niayesh Afshordi and Stephen Ro 



log, (p) [gr/cm 3 ] NO COOLING 



lo gioC°) [g r / crr|3 ] cooling 



= 0.8Gyr 




= 0.8Gyr 





.6Gyr 




t=1.6Gyr 




>- 200 



t = 2.4Gyr 





3.2Gyr 



t = 3.2Gyr 




: 










i i i 




: 













100 


200 


300 


400 500 
X(kpc) 


600 700 800 


900 





100 200 


300 400 500 
X(kpc) 


600 700 800 


900 




-28.0 




7.5 


-27.0 


-26.5 


-26.0 -25.5 


-25.0 


-28.0 


-27.5 


-27.0 -263 


-26.0 -253 


-25.0 



Figure 2. Slices of gas density (logarithmic) in the XY plane through the middle of the simulation box for the non-radiative (left) and 
dissipative (right) simulations at t = 0.8, 1.6,2.4 and 3.2 Gyr. 



does not cool further. We remark that this feedback does 
not affect the regions of interest, namely outside ~ 0.7r v 
(the stripping radius) and that we do not intend to give a 
realistic feedback model but simply to attenuate the central 
overcooling. We simulate a single gas species with a metal- 
licity of \og\o(Z / Zq) — —1.5 in both, the satellite and the 
ICM. Although in a realistic setting the gas is likely to have 
different metallicities according to its position relative to 
the center of the satellite, and to the center of the low-mass 
group, we choose this value for simplicity and because it fur- 
ther help us reducing the central collapse of the corona. We 
note that this value is roughly of the order of what is found 
in simul ations of MW-size haloes near the virial radius (see 
Fig. 9 of Ivan de Voort fc Schavell2012l ). 



4 RESULTS 



Fig. ((2j) shows slices of the XY plane through the middle of 
the simulation box for the density of the gas for the non- 
radiative (left panels) and dissipative (right panels) cases at 
four output times t = 0.8, 1.6, 2.4 and 3.2 Gyr from top to 
bottom, respectively. The effect of RPS and the development 
of KH instabilities are clearly seen. In the top panels, we see 
how the shock has propagated through the satellite and is 
passing through its dense centre, whereas in the outskirts, it 
has traversed most of the satellite. By t ~ 2.4 Gyr, the shock 
has crossed the satellite completely and the gas removed by 
RPS is already forming clear trailing wakes. This roughly 
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agrees with our estimate of the RPS time scale to be of the 
order of 1.8 Gyr (see Table [2] and Appendix |A"| . 

In the dissipative simulation, the central collapse is al- 
ready apparent after t ~ 1 Gyr with the formation of a dense 
and cold core. Although the artificial suppression of cooling 
prevents the central collapse from continuing further, the 
contraction of the corona is large enough to significantly re- 
duce the amount of stripped gas. However, the presence of 
the trailing wakes and the effect of turbulence are still evi- 
dent. The linear, dense and cold feature to the right of the 
center of the satellite is caused by the asymmetry of the 
problem and the inclusion of a spherically symmetric cool- 
ing suppression. Although the shape of the central "inner" 
wake is affected by this feature, there is only a minor impact 
on the extended "external" wakes that participate in most 
of the mixing with the ICM. 

4.1 Turbulence and gas mixing 

After t ~ 1 Gyr, the stripped gas starts mixing efficiently 
with the ICM through the turbulence generated by KH in- 
stabilities. From an initial look at the density and tempera- 
ture of the gas at t — 3.2 Gyr, it seems that the gas mixture 
has twice the density and nearly half the temperature of the 
ambient ICM. To give a first estimate of the amount of mass 
that makes the bulk of the mixed phase, we show in Fig. ([3]) 
a mass-weighted temperature histogram of all the mass con- 
tained within a box that extends above and below the slice 
shown in Fig. ([2]) with a thickness of ~ 490 kpc. In this way, 
the volume covered by this box encloses all the mass that 
is being mixed, the satellite, and of course a fraction of the 
ICM that is not involved in the mixing. The black histogram 
is for t = Gyr, while the blue and red are for t = 3.2 Gyr 
for the simulations with and without cooling, respectively. 
The initial gas distribution within the satellite is clearly seen 
around 10 6 K as a "cold" phase, as well as the singled- value 
"hot" phase in the ICM to the right. At the end of the non- 
radiative simulation, the gas in the satellite has been shock- 
heated while the ICM has been adiabatically compressed 
and expanded (see below). For the dissipative simulation, 
the centre of the satellite has partially collapsed, forming a 
cold core with a temperature given by the imposed floor of 
10 s K. The rest of the gas however seems to be distributed 
in a similar way in both simulations; in particular, the for- 
mation of a "warm" phase at T ~ 2 x 10 6 K is apparent in 
the figure for both cases. We show below that this phase is 
actually formed by the mixture of the two phases and not 
just by the gas stripped from the satellite. 

The upper panels of Fig. Q give an idea of the degree 
of mixing between the fluids at the end of the non-radiative 
(left panel) and dissipative (right panel) simulations. They 
show the distribution of tracer particles that were set at 
t = to sample the density distribution in the satellite (red) 
and in the ICM (black). In the figures, a projection of this 
distribution at t = 3.2 Gyr is shown for a slab with a thick- 
ness of 50kpcB The objective of this figure is to show that 
the two fluids originally separated are thoroughly mixed by 

6 The tracer particles sample the gas distribution with lower res- 
olution than the mesh cells, particularly for the ICM where the 
initial inter-particle separation is ~ 15 kpc. 
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Figure 3. Mass-weighted histogram of the temperature within a 
reduced volume of the simulation box, see text for details. The 
black histogram is for the initial simulation time, whereas the red 
and blue histograms are for the final simulation time for the cases 
without and with radiative cooling, respectively. 



the turbulent motions in the plumes. The degree of mixing 
is similar in both simulations. 

The bottom panels of Fig. Q are a slice through the 
middle of the simulation box showing the cooling time, con- 
structed using the temperature and density of the gas at t = 
3.2 Gyr and assuming a single metallicity, \ogio(Z/Zo) = 
—0.5, representative of the intragroup gas. In reality, the 
cooling time of the gas mixture would depend on the mass- 
weighted metallicities of the species being mixed. Since these 
are uncertain, and for simplicity, we have chosen a constant 
metallicity, but of course, the absolute value of the cooling 
times reported in this work are metallicity-dependent. We 
see that in the trailing wakes where the gas mixes more effi- 
ciently, the cooling times have decreased by up to an order 
of magnitude compared to the original cooling time in the 
ICM. 

Fig. (0 shows the position of the tracer particles (in red 
and black for the satellite and ICM, respectively), within the 
slab shown in the upper panel of Fig. Q (i.e., within a box 
of volume 942 x 487 x 50 kpc 3 ) in a temperature-density 
diagram at different simulation times: t = 0.1,0.8,2.4 and 
3.2 Gyr clockwise starting from the upper left of each set of 
panels. The set to the left (right) corresponds to the sim- 
ulation without (with) radiative cooling. The dotted lines 
are lines of constant cooling time from 100 Gyr to 0.1 Gyr 
going from left to right. The thick black solid line shows 
the locus of points where the cooling time equals the local 
free-fall time (tg = i C ooi). Finally, the dashed line is a refer- 
ence marking processes that conserve entropy (defined here 
as S ocT/p 2/3 ). 

Since the number of tracer particles in the satellite is 
high (particularly in the center of the satellite), we have 
randomly selected a smaller sample to plot in Fig. JSJ, ex- 
cept in the region around the wakes for the last two time 
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Figure 4. Upper panels: Distribution of tracer particles after t = 3.2 Gyr for a slab through the middle of the simulation box with a 
thickness of 50 kpc; red and black for the satellite and ICM, respectively. Bottom panels: Slice in the XY plane through the middle of 
the simulation box at t = 3.2 Gyr showing the values of the cooling time (logarithmic) assuming \og-K)(Z / Z q) = —0.5. The panels on 
the left (right) correspond to the simulation without (with) radiative cooling. 



outputs, bottom panels, where the tracers are shown with 
bigger symbols. We have verified that this random subsam- 
ple traces the distribution of the full sample and remark 
that this was done simply for plotting purposes. This region 
around the wakes is the one corresponding to X > 500 kpc 
in Fig. (see section g^ . 

Initially, the ICM gas represents a single point in this di- 
agram (p = 10" 28 g cm" 3 , T = 3.6 x 10 6 K) and the gas in the 
satellite is distributed according to the equilibrium configu- 
ration. In the first 100 Myr, the gas in the ICM immediately 
surrounding the satellite compresses upstream and expands 
downstream, conserving its entropy. The gas on the satellite 
is compressed and shocked upstream; the gas loosely bound 
to the satellite begins to be removed. As the shock traverses 
the satellite, it heats and compresses the gas within it, re- 
moving material from the outskirts while gravity tends to 
restore the original gas distribution; the distribution within 
the inner regions will remain distorted though, since the 
shock causes an irreversible change in the gas. 

In the dissipative simulation, the gas that remains 
within the satellite is additionally subjected to the central 
collapse, which is suppressed in the innermost regions by 
the feedback-like mechanism we have introduced. This cre- 
ates a gas distribution that is highly dispersed in the phase 
diagram (the coldest and densest material is not explicitly 
shown in Fig. [S] since we concentrate on the regions where 
mixing happens). On the other hand, the behaviour of the 
gas in the outskirts of the satellite, and outside in the ICM, 
is qualitatively similar to the one in the non-radiative simu- 
lation. In particular, most of the ICM evolves along a path of 



constant entropy (expanding and contracting driven by the 
relative motions of the fluids, and the pressure support and 
gravity of the satellite). At t ~ 1 Gyr, a fraction of the ICM 
is clearly mixing with the stripped gas increasing its density 
and decreasing its temperature. By t — 3.2 Gyr, a mixed 
phase of gas is apparent in the figure with T ~ 2 x 10 6 K 
and p ~ 2 x 10~ 28 g cm -3 . This confirms the initial visual 
impression given by the bottom panels of Fig. @, and also 
proves that the "warm" phase shown in Fig. © is a mixed 
phase formed by gas from the satellite and the ICM. 

By looking at the lines of constant cooling time, we 
see that the mixed phase reaches values of t coo \ ~ 5 Gyr, 
considerably lower than the initial cooling time of ~ 40 Gyr. 
We also see that by the end of the simulation, the mixed 
phase is approaching the locus of points ts = icooi that 
separates the gas that can cool efficiently from the one that 
cannot. 



4.2 Properties of the wakes 

We define the wakes as the trailing structures that are clearly 
distinguished in Fig. ((2| at X > 500 kpc. By doing so, we 
exclude the material that is not bound to the satellite but 
is still attached to it. This gas will be stripped in the future 
and added to the trailing wakes, but since this gas is clearly 
not participating in the mixing we omit it from the analy- 
sis. We note that although this definition is arbitrary, the 
properties of the wakes at a given time serve as an approxi- 
mate reference for the properties at other times. With such 
definition, we find that at t = 3.2 Gyr roughly 50%, 55% of 
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Figure 5. Temperature-density diagram for the tracer particles contained within the volume that appears in the upper panel of Fig. [4] 
Only a random subsample of the full tracer population within this region is plotted (except in the bottom panels for the region around 
the wakes, X > 500 kpc in Fig. [2] where the tracers are shown with bigger symbols. The red (black) points are for the satellite 
(ICM). The simulation times shown are t = 0.1,0.8,2.4 and 3.2 Gyr. The dotted lines are lines of constant cooling time (assuming 
logio(^/^0) = —0.5) and the black solid line shows the locus of points where % = t coo i. The dashed lines mark an example of a path 
with constant entropy. The set of panels to the left (right) corresponds to the simulation without (with) radiative cooling. 



the mass that is no longer bound to the satellite is in the 
trailing wakes in the simulations with and without radiative 
cooling, respectively. In the former the amount of stripped 
gas is ~ 7 x 10 9 Mq, a factor of 3 smaller than in the latter. 



The gas in the ICM mixes with the gas from the satel- 
lite and becomes part of the trailing wakes. Although this 
mixing occurs in a wide region, we can arbitrarily define 
the "warm phase" as the region in the wakes where the 
gas (both from the ICM and from the satellite) has an 
entropy which is at least 50% lower than the initial en- 
tropy of the ICM. Gas with higher entropy has either not 
been thoroughly mixed yet and/or contains the highest en- 
tropy gas from the satellite. Using this definition, we can 
narrow down the median values of the temperature and 
density of the gas that is more likely to cool afterwards 
and estimate its total mass after 3.2 Gyr of the evolu- 
tion for the simulation without (with) radiative cooling: 
<T mix ) ~ 2(2.1) x 10 6 K, (p mix ) = 1.8(1.8) x lCT 28 g cm" 3 , 
and M mix = 7.6(1.4) x 1O 9 M (~ 20(28)% coming from the 
ICM). These values imply i coo i ~ 7(8)Gyr, a factor of ~ 6(5) 
less than that of the ambient ICM. In the non-radiative 
simulation, most of the gas in the wakes that came from 
the satellite belongs to the warm phase (59%), while in the 
dissipative simulation is a minority (27%). This is simply 
explained by the smaller amount of stripped gas that par- 
ticipates in the mixing in the dissipative simulation, which 
also explains why the warm phase has a longer cooling time 
in this case. 



5 DISCUSSION 

The relatively weak impact of radiative cooling in the for- 
mation and initial properties of the wakes (with the main 
difference being a reduction of the stripped gas) can be un- 
derstood by noting the following. The gas in the satellite 
that participates in the mixing is given, at the most, by the 
gas within a shell external to the stripping radius. Even if 
all the gas within this radius were to collapse immediately 
due to radiative losses, the outer shell would collapse in a 
free-fall time, which is > 1 Gyr; this contraction would hap- 
pen adiabatically since the gas in the shell has still quite a 
long cooling time. An adiabatic process preserves the cool- 
ing time at the temperatures and densities of the collapsing 
shell (note how the lines of constant entropy are roughly 
parallel to the lines of constant cooling time in Fig. [5)l . This 
is because the cooling function goes as A coo i oc T~ a , with 
a ~ 0.5 for logio(Z/Z@) ~ —0.5, in the temperature range 
between 10 5 K and ~ 5 x 10 6 K. Thus, the gas in the shell, 
that is eventually stripped, has enough time to mix with the 
ICM forming a wake in a qualitatively similar way as in the 
non-radiative simulation. 

In detail, the properties of the wakes depend on the 
evolutionary stage of the interaction and, of course, on 
the metallicity, temperature and density of the gas being 
stripped. The latter depend on the assumed radial profiles 
and have important consequences in our results. Although 
we cannot check directly from observations if these profiles 
are adequate, we found that the the properties of the gas 
near the outskirts of the corona are ro ughly the same as 
the ones reported in ICrain et all (|2010h (based on a large 
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sample of MW-size galaxies extracted from a cosmological 
simulation). The mass of stripped gas not only depends on 
the balance of radiative cooling and central feedback (which 
determines how fast the central collapse relative to the strip- 
ping time scale is) but also in the uncertain amount of hot 
gas that is in the corona. In this work we have used the 
largest possi ble value set by the un i versa l bary on fraction. 
However, e.g. lAnderson fc Bregmanl (|201p| . l201lh found that 
the actual amount of mass in the hot coronae is consider- 
ably less than naively expected from subtracting the uni- 
versal baryon fraction from the amount of baryons that is 
locked in stars and cold gas. A lower gas fraction would 
lower the mass of stripped gas. Because of all this, our pre- 
dictions cannot be made in terms of the absolute values of 
the properties of the gas that mixes, but rather in terms of 
the mixed gas fraction relative to the stripped gas, and its 
cooling efficiency relative to the temperature and density of 
the removed gas. 

Our non-radiative simulation indicates that the analytic 
model we introduced in Section [2] gives reasonable estimates. 
In particular, once we apply the model to the conditions in 
the simulation, we find that it is able to reproduce some of 
the global properties of the wakes (see below) once we adjust 
its parameters and add two additional components to it: i) 
The stripped gas is shock-heated to temperatures in excess 
of the virial temperature of the satellite before it mixes with 
the ICM. We can account for this effect by calculating the 
Rankine-Hugoniot (RH) temperature jump conditions of the 
shock (leading to a shock-heated temperature ~ 1.7 x 10 6 K); 
ii) The gas in the wakes is compressed by a factor <5 m i x rel- 
ative to the density of the ICM. With these considerations 
we find that by taking: a rem ~ 1, ft ~ 3, and <5 m i x ~ 2 we 
are able to reproduce the amount of stripped gas, the me- 
dian temperature in the mixed phase and its cooling time. 
Although comparing the analytical model to the dissipative 
simulation is in principle more complicated due to the con- 
traction of the corona, it is possible to use the current model 
without modifications if we assume that the contraction oc- 
curs rapidly relative to the gas removal processes. Then, the 
results of the dissipative simulation can be matched by set- 
ting Qfrcm = 1/3 and /3 — 1.5. The first of these parameters 
is reduced, relative to the non-radiative case, to accommo- 
date a smaller amount of stripped gas, while the second is 
reduced since the corona is more compact. 

Perhaps the most important limitation of the model 
arises from the complex nature of the mixing process which 
occurs in the turbulent eddies in layers of different size. This 
implies that defining an overall effective value for /3 is unreal- 
istic since there is a broad distribution of cross sectional ar- 
eas. In the simulations, the mixing happens along the annu- 
lus formed by the wake, which has a cross sectional area that 
broadens if the relative velocity of the satellite decreases. 
Although the gravitational focusing effect is stronger in this 
case, the gas in the turbulent eddies is confined to a re- 
gion which is more extended than the case of larger velocity 
since the ram pressure of the surrounding ICM is smaller 
(we have corroborated this by running an identical simu- 
lation but with i; S at = 150km/s). This suggests that j3 in- 
creases (decreases) for lower (higher) velocities, contrary to 
expectations. 

Despite its simplicity, the analytic model could be used 
for instance in semi-analytic models (SAMs) of galaxy for- 



mation (for a review see e.g. lBaughll200fj| ) to account for 
the removal processes of the coronal gas, which are typically 
severely simplified by assuming that this gas is removed in- 
stantaneously having no further importance. 

There are specific conditions one can think of that 
would enhance the mechanism proposed in this work. For 
instance, if the infalling satellite has lower mass and hence 
a colder corona, the temperature in the mixture would be 
lower than the case we analysed, and gas removal would be 
more efficient. A lower infall velocity would create a weaker 
shock and thus a reduction on the temperature of the shock- 
heated gas as it is stripped from the satellite On the other 
hand, the strength and time scale of gas removal decrease 
making the mixing process less efficient. 

Although more speculative, we also expect that the 
mechanism would be more important as the satellite spi- 
rals towards the centre of the dense stratified medium of 
the host. We can crudely estimate what would happen in 
these circumstances using our analytic model. Although the 
process is of course gradual, let's assume that by the time 
the satellite is at a distance of 0.5r Vi host from the centre of 
the host, it has lost an amount of dark matter and gas as pre- 
dicted by the analytic estimates at the virial radius, i.e., the 
satellite arrives at this orbital distance with a mass profile 
"clipped" at a tidal stripping radius r t = 0.65r v (see Table 
|2J. If the specific angular momentum of the satellite was 
conserved along its orbit, then its velocity would be twice 
the velocity at the virial radius. Dynamical friction how- 
ever, decelerates the satellite; for simplicity we assume it to 
be ~ l.Susacl The temperature in the corona of the satellite 
is that left by the remnant of the initial shock, which is given 
by the RH jump conditions of the virial shock. The temper- 
ature and density of the ICM at this position are 3.5 x 10 B K 
and 6 x 10 _28 g cm -3 , respectively. Under these conditions, 
we applied our analytic method and find that the remaining 
corona would be removed in a time scale roughly given by 
the sound crossing time < 1 Gyr, and the gas in the mixture 
would have a cooling time of only 1.5 Gyr (a factor of ~ 4 
less than the ambient ICM). 

It is worth mentioning that recent hydrodynamical 
simulatio ns of galaxy fo rmation performed with the code 
AREPO (ISprin gcl 2010a, a hybrid method combining some 
of the advantages of Lagrangian and Eulerian codes and 
avoiding their most important weaknesses) have reported 
significantly higher star formation rates for central galax- 
ies in massive haloes than the ones found, under the same 
conditions, in an SPH formulation l|Keres et al.ll201ll ). The 
authors advocate that the difference is partly due to a more 



The typical satellite infall velocities are of the order of the virial 
velocity o f the group, lo wer infall velocities are nevertheless pos- 
sible (e.g. Wetzel 201l|). For example, if u sa t = 150 km/s, then 
the temperature ratio of the post- to pre-shock gas according to 
the RH temperature jump conditions is reduced from 2.8, in the 
case when Vaat = 315 km/s, to 1.6. 

8 We ran a dark-matter only simulation of a MW-size halo merg- 
ing with a 10 13 Mq group. The initial infall velocity of the satel- 
lite, with its radial an d tangential co mponents, was set to the av- 
erage value as given in Wetzel (2011). We found that the satellite 
reaches an orbital distance of 0.5r v host with a speed of ~ 1.6i) sa t 
with approximately 70% of its original dark matter mass. This 
simulation confirms the orbital evolution we have assumed here. 
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efficient mixing in AREPO, caused by a l arger removat of fow 
entro py materiai f rom infalling satellites (|Vogelsberger et al.l 
l201lf ). Fig. 4 of iKeres et all (|201lh shows that for a ~ 
1O 13 M0 halo at the present day, this difference is a factor 
of ~ 2: AREPO predicts approximately 3 x IO^Mq of ad- 
ditional stars. It is unlikely that the majority of this extra 
mass can be accounted for by the coronae mixing studied in 
this paper, however it might be that a significant portion of 
the additional star formation material comes from a similar 
mixing process between the ICM and the cold gaseous disc 
of the merging galaxy. We believe it would be interesting to 
explore the relative contribution of both of these processes. 



6 CONCLUSIONS 

The presence of an extended corona of hot gas within clus- 
ters, groups and massive galaxies is a generic prediction of 
the current theory of galaxy formation. Since such coronae 
are a large gas reservoir that can potentially cool down and 
form stars, it is of great relevance to study their evolution 
and in particular their interaction during galaxy mergers 
which are ubiquitous in the CDM paradigm. 

Most studies of galaxy mergers have concentrated on 
the interaction of the cold components (stars and gaseous 
discs), relegating the hot coronae to a secondary role es- 
sentially governed by ram pressure stripping. In the case 
of MW-size galaxies merging into groups, the high tem- 
perature and low density of the ambient ICM at the out- 
skirts of the groups are commonly assumed to imply that 
the stripped gas stays irrelevant for subsequent cooling and 
star formation. This assumption, however, ignores the fact 
that Kelvin-Helmholtz instabilities produced by the velocity 
shear also remove coronal gas and efficiently mix it with the 
ICM through turbulent eddies. Since the former is colder 
than the latter, the mixture could have a shorter cooling 
time than the ambient intragroup gas. These instabilities 
are typically suppressed in the majority of galaxy formation 
simulations that have been done to date since most of them 
were carried out using the SPH formulation l|Agertz et al.l 
I2007I ). 

In this paper, we use analytic methods and hydrody- 
namical simulations (with and without radiative cooling) to 
study the gas removal processes that act during the inter- 
action a of a spherical MW-size corona moving through a 
homogeneous ICM wind representative of the coronal gas at 
the virial radius of a low-mass group. The main conclusions 
of our work are the following: 

• Our simulations show that a relevant fraction of the 
removed coronal gas mixes with the ICM in trailing wakes 
creating a "warm phase" (i.e. with at least 50% less entropy 
than the ambient ICM) with a cooling time which is almost 
an order of magnitude lower than the cooling time of the 
ambient intragroup gas. This reduction is a result of a lower 
temperature in the gas mixture driven by the colder coronal 
gas from the satellite, as anticipated, but also due to an 
enhanced density caused by compression. 

• We find that the gas mass in this warm phase is between 
40-76% of the mass stripped from the satellite and deposited 
in the wakes (M war m ~ 1.4 — 7.6 x 1O 9 M0 for the parameters 
we have used in our simulations). This gas is mostly com- 



ing from the satellite, but a non-negligible fraction, 20-28%, 
comes from the ICM. 

• The analytic model described in Section [2] is able to 
roughly estimate the amount of removed mass and the time 
scale of removal due to the two main acting mechanisms: 
ram pressure stripping and KH instabilities. Moreover, once 
properly calibrated, it can be used to estimate the most 
likely temperature and cooling time of the gas mixture af- 
ter the gas settles in the turbulent wakes. We believe that 
this model can be incorporated into semi-analytic models of 
galaxy formation to account for these processes, which are 
typically ignored in current implementations. 

We speculate that the warm phase created through the 
mechanism proposed in this work has the potential to be- 
come a mode of star formation in galaxy groups (through 
in situ star formation and/or through a cold stream flowing 
to the central galaxy of the host). In fact, the possibility of 
stars forming in turbulent wakes from the condensation of 
gas stripped from cold discs seems to be a likely hypothesis 
in some observed cases (e.g. iHester et al.ll2010l ). Although 
in such cases the mechanism at work is the same we have 
studied here, star formation is far more likely in those since 
the gas being removed is much colder than the ambient gas. 
We are thus cautious in speculating about the fate of the 
removed coronal gas. In particular, it is not clear whether 
the gas mixture can survive long enough without evaporat- 
ing th rough conduction with the ICM fe.g. lNipoti fc Binnevl 
1 20041 ). This would depend on the specific details of conduc- 
tion in the surfaces of the turbulent wake, which are difficult 
to estimate. It is therefore crucial to address these questions 
using more detailed simulations. We believe that the results 
of this work motivate further analyses to fully quantify the 
effect of the mixing of the coronal gas from satellites with 
the intragroup gas. 
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APPENDIX A: RPS TIME SCALE 

The properties of the initial shock generated by the collid- 
ing coronae and the values of the post-shock thermodynamic 
variables can be estimated by solving a one-dimensional 
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shock tube problem. Once the shock develops, let region 1 to 
the left be the unperturbed ICM moving to the right with 
a relative velocity w sa t and with properties (p, P, u, c s )i = 
(Phost, fhost, «sat, c s ,host), and region 4 to the right be the 
unperturbed corona of the satellite at rest having proper- 
ties: (p,P,u,c B ) 4 = (p sa t(r v ),Psat(r v ),0, c SiSat (r v )). Regions 
1 to 4 separate the three characteristics of the problem. The 
forward shock front divides regions 3 and 4; thus, region 3 
contains the post-shock material. The solution to this prob- 
lem is found by solving the following transcendental equa- 
tion: 



Pi 

P4 
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i + 
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where (, — P3/P4 and A = (7 + l)/2 r y, and we have assumed 
that the corona of the satellite and the ICM are in pressure 
equilibrium at first: P1/P4 = 1. The velocity of the shock, 
and the velocity of the post-shock fluid are then given by: 
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Once the shock is generated it will propagate towards the 
centre of the satellite. The density stratification of the lat- 
ter will then change the speed of the shock and modify the 
post-shock properties of the newly encountered material. To 
roughly estimate the propagation of the shock, we follow 
the description of iKogure fc OsakH (|l962h that reduces the 
evolution of the post- and pre-shock pressure ratio £ to the 
following ordinary differential equation (derived from a one- 
dimensional problem neglecting the effects of gravity and 
radiative losses, see their Eq. 3.1): 



_ (1 + Q(l + AQ(2 + (3 - A)Q rilnP 
ds' ~ 4 + 2(2 + A)£ + A(3- A)£ 2 ds' 
^ g(l + C)(l + AC) dlnp 



4 + 2(2 + A)C + A(3-A)0 ds' 



(A3) 



where s' is the Lagrange coordinate of the shock front. In 
reference to the variable s = r/r v from the centre of the 
satellite in Eqs. {1} and p|. s' — 1 — s. Thus, solving 
Eq. (|A1|) assuming pressure equilibrium in the boundary 
we get s ' — 0) = 9.0(12.9) for the galaxy-group (galaxy- 
cluster) interaction. With this boundary condition, we can 
solve Eq. (|A3|) to find £(s') and use it in the first of Eqs. (|A2[) 
to get the shock speed as a function of s' , which can then 
be substituted in Eq. (JS| to find the time scale for RPS. 
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